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Abstract. 

We present preliminary results of our implementation of a novel electrophoresis separation tech- 
nique: binary oscillatory cross flow electrophoresis (BOCE). The technique utilizes the interaction 
of two driving forces, an oscillatory electric field and an oscillatory shear flow, to create an active 
binary filter for the separation of charged species. Analytical and numerical studies have indicated 
that this technique is capable of separating proteins with electrophoretic mobilities differing by 
less than 10%. With an experimental device containing a separation chamber 20cm long, 5cm 
wide, and 1mm thick, an order of magnitude increase in throughput over commercially available 
electrophoresis devices is theoretically possible. 

Introduction. 

Electrophoresis has long been shown as an effective process for the separation of particles and 
biological macromolecules based on small differences in electrophoretic mobility. Since Tisilius [1] 
discovered electrophoresis in 1937, the technique has become a cornerstone for biochemical analysis 
and has branched into several different subfields. Despite its attention, several key problems have 
prevented electrophoresis techniques from scaling up to commercial processes. Among these are 
Joule heating and crescent dispersion caused by electroosmosis, parabolic velocity profiles of the 
carrier electrolyte, and electrohydrodynamic effects [2], [3]. 

The present research project, binary oscillatory crossflow electrophoresis, has the potential for 
significant improvements in throughput and resolution over conventional CFFE techniques. This 
technique relies on the interaction of an oscillatory electric field and an oscillatory shear flow, shown 
in Figure 1. 
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Figure 1: Direction of motions in a binary oscillatory crossflow electrophoresis cell. 


By appropriate selection of frequency and phase of oscillation of the two interacting driving forces, 
an effective filter can, in principle, be created that will allow proteins either higher or lower than 
a chosen electrophoretic mobility to pass through the device. The technique could be used for 
both analytical and preparative work by programming the filter setpoint with time and collecting 
proteins that pass through the device at each setpoint. For continuous isolation of a particular 
fraction from a complex mixture feed, a network of four of these devices could be used. 

Analytical Model 

The driving forces in the binary separator are an oscillatory convective flow and an oscillatory 
electric field. The solute molecules will have some net motion that is a complex function of the two 
driving forces, the electrophoretic mobility, and the diffusion of the species. The cross flow may be 
written as 


u x = Au max y - - |) f(ut) 
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where u x is the fluid velocity in the x direction, u max is the velocity of the center line, y is position 
in the gap, b is the gap thickness, u is the frequency of oscillation, t is time, and / is some periodic 
function of ut. The motion of a particular solute molecule across the channel will depend on /(u/£) 
and the position of the molecule in the gap. If we neglect diffusion, we may write the motion of 
the solute across the gap due to an oscillatory electric field as 
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where /.i is the electrophoretic mobility, E 0 is the characteristic amplitude of the electric field, g(ut) 
is the oscillatory motion of the solute, and h(ut) is the oscillation of the electric field. Integrating 
this equation, we can solve for the position across the gap of the solute molecule. 

| = ag(ut) (3) 

where a = is the dimensionless electrophoretic mobility and g{wt) is bounded between 0 and 
2. Substituting the position of the solute molecule into equation 1 and taking the time average, we 
will have an equation that describes the motion of a solute molecule in the absence of molecular 
diffusion: 

( ^f) = 4UmaxQ (<0M)/M)> ~ (flM)VM))) ( 4 ) 

where (•) is the time average over one period. Equation (4) is limited to dimensionless mobilities 
a < 0.5. The maximum value of a corresponds to a solute molecule that just reaches the opposite 
wall during one period of oscillation. With equation (4) we may evaluate the solute velocity across 
the channel for combinations of arbitrary choices for f(u>t) and g{u)t). Chandhok and Leighton [4] 
chose g{u)t) = 1 — cos(uf) and showed that the only Fourier modes of a general periodic function 
f{vt) that contributed to a migration velocity were a steady flow, cos(u;£), and cos(2u>£). While this 
analytic solution is very useful for qualitatively determining the performance of such a separations 
device, the solution is seriously limited by neglecting the effects of diffusion. This is particularly 
true since the interaction of diffusion with an oscillatory shear flow has been shown to increase the 
dispersion in the direction of motion by several orders of magnitude [5], [6]. Furthermore, bounding 
the maximum value of a significantly limits the possible operating parameters of the system. 

In order to better characterize the motion of the solute, we numerically simulated the solute 
motion in the channel including the effects of diffusion. To perform the numerical simulation, we 
used a molecular dynamics approach. At each time step, we added a random walk of length \/2AA t 
to the electrophoretic motion of a solute molecule and used the resulting position to calculate the 
displacement in the direction of the fluid motion. The simulation was carried out for 1000 molecules 
and the average position and variance of the solute as a function of time were recorded. From the 
position data we determined the time averaged solute velocity and the effective diffusion coefficient 
in the direction of motion of the solute. We ran the simulation until the velocity reached a steady 
value, usually after only a few periods of oscillation. 

For our simulation we have chosen the gap width b = 0.1 cm, and by arbitrarily choosing the 
dimensional electrophoretic mobility of a particular species and the oscillation frequency u, we 
fix the electric field amplitude for a desired value of a. Consequently, our choice of the above 
parameters fixes the dimensionless diffusivity. We simulated the solute motion for an electric field 
of amplitude h(t) = sin (u>t) - 0.05 coupled with both f(u>t) = cos [u>t) and f{u>t) = cos(2u;£). The 
small steady component to the electric field strength of 0.05 is important to refocus the solute at 
the lower wall during each period. The velocities for the solute as a function of a for f(t) equal to 
cos(utf), cos(2u;£), and cos{u;t) - cos(2u>f) coupled with h(t) — sin (ut) - 0.05 were simulated and 
are shown in Figure 2. 
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Note that f(t) equal to cos(wt) and cos(2wf) both lead to negative velocities, with the maxi- 
mum amplitude occurring at different mobilities. For a binary separation to be possible, the time 
averaged velocity for the solute molecules must change sign at some critical value of the elec- 
trophoretic mobility. Because the fluid motion does not directly influence the migration across the 
gap, the effect of each mode of the imposed solvent velocity on the time averaged solute velocity 
is linearly additive. As a result, we chose to subtract the two modes. The velocity resulting from 
f(u>t) = cos(u;t) - cos(2 ut) coupled with h(t) = sin (u>t) - 0.05 meets the constraints for a binary 
separation. In addition, j s large near the mobility with zero net velocity, a condition for 

sharp separation. 

To determine the characteristic throughput and selectivity of the device, we consider a cell 
connecting two reservoirs of fluid. In one reservoir, we impose some concentration cai, in the other 
reservoir we impose a concentration of zero by flushing it with buffer. The cell is of length L, width 
W, and thickness b with operating parameters as listed above. The cross flow will oscillate with 
f(ut) = cos(u;t) - cos(2 lot) and u max = Axcj where Ax is the characteristic amplitude of fluid 
oscillation. 

To find the flux though a cell we start with the governing equation 
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where k e jj is the dispersion coefficient determined numerically, x is the direction of flow, ca is the 
concentration of species A, and va is the velocity of species A depicted in Figure 2. Equation (5) 
is the result of substituting the flux into the mass balance. By integrating twice and imposing the 
boundary conditions of c^(x* = 0) = 1 and c^(x* = Lj Ax) = 0, we can solve for the dimensionless 
concentration profile. Substituting the dimensionless concentration profile into the flux equation, 
we find the the dimensionless flux is 


N'a = V A ' 


1 — exp 



( 6 ) 


Figure 3 shows a plot of as a function of a for a value of ^ | where u \ is as shown 

in Figure 2 and k*jj is determined numerically. The curve labeled high pass represents the flux 



Figure 2: Plot of dimensionless velocity versus dimensionless electrophoretic mobility including the 
effects of cross-streamline diffusion. 
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when f(uit) = cos(u >t) - cos(2wf), and the curve labeled low pass represents the flux when f{ut) = 
cos(2u;t) - cos(wt). 



Figure 3: Plot of Dimensionless Flux versus dimensionless electrophoretic mobility. 


Clearly the flux changes very dramatically near a — 0.41 for both high and low pass operation. 
For a 10% change in a about this point, the flux increases two orders of magnitude. 

Experimental. 

To verify the analytical and numerical results of the BOCE technique, an experimental system 
including a single binary separation cell, two reservoirs and an electrolyte bath was constructed. 
The cell is similar to that of Giddings [7], and allows continuous flow of solute in a buffer solution 
through a separation chamber 20.0 cm long, 5.0 cm wide, and 0.1 cm thick. Care was taken to 
address problems inherent in electrical flow cells of this type; mainly gas generation from electrode 
reaction, Joule Heat and membrane flexure. 

The BOCE cell is approximately 24.0 cm long, 7.0 cm wide, and 2.0 cm thick. The central 
separation chamber of the cell is machined from a plastic sheet that is divided into two halves to 
accommodate feed and exit channels. Fluid is introduced at either end of the chamber by means of 
a stainless steel syringe connected to polypropylene fittings extending through the electrolyte bath 
and connecting to the reservoirs. The separation chamber is sandwiched between two cellulose 
membranes, supported by the electrode screens. The layers are clamped between two plexiglass 
blocks. In the center of each block a rectangular opening approximately the size of the separation 
chamber was machined. This design allows the central chamber to be sealed with the membranes 
and screen support using bolts to provide uniform pressure around the perimeter. 

Giddings [7] found in experiments using a similar device that cellulose membranes are prone to 
flexing even upon small pressure fluctuations. The requirement of fixed volume in the separation 
chamber thus necessitates the use of a support. Aside from chemical compatibility with the buffer 
solution and mechanical strength, the electrical resistance was an important consideration. It was 
therefore decided that the electrode screens could be used as supports, thus serving dual purposes. 
This reduced the electrode separation to under 1.1 mm. 

Gas bubbles and Joule heat are removed by submerging the entire cell in an electrolyte bath 
open to atmosphere. The bath accommodates up to 6 liters of electrolyte fluid and copper cooling 
coils. The cell is positioned in the bath with the 7.0 cm side forming the vertical axis, and gas 
bubbles generated from electrode reaction are removed from the electrodes through the machined 
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openings in the clamping block halves. 

The two oscillatory driving forces were provided by a KEPCO BOP100-4M bipolar power supply 
capable of constant current operation and a specially designed syringe pump. The electric field 
was controlled using LabVIEW software and a 100 MHz Pentium personal computer. Bulk fluid 
conductivity data was acquired into LabVIEW using a Cole-Parmer conductivity meter and the 
amplitude of the current waveform was simultaneously adjusted to maintain the desired electric 
field. The flow waveforms were delivered to the syringe pump using a Galil DMC-1500 motion 
controller and software. Position data from the syringe pump was fed back into LabVIEW allowing 
phase locking of the two driving forces to within a few milliseconds. Temperature was monitored 
using the Cole-Parmer conductivity meter and protein concentration was analyzed using a Varian 
UV-VIS spectrometer. 

Results 

The focus of preliminary experiments was to verify qualitatively that the interaction of a purely 
oscillatory shear flow with an oscillatory electric field could indeed lead to the net convection of a 
single charged protein species under conditions of uniform concentration throughout the system. 
Furthermore, we sought to demonstrate that by changing the dimensionless mobility, a, by changing 
the amplitude of the electric field, we could change the sign and magnitude of the velocity of the 
protein. To this end preliminary experiments were conducted using Bovine Hemoglobin (BHb) 
dissolved in a sodium acetate/acetic acid buffer of pH 4.5 and conductivity of approximately 
375.0 Based on results from Douglas et. al. [8], the electrophoretic mobility of BHb was 
approximately 17.4e _5 ^j. The choices of waveforms included h(wt) = sin(u>£) - 0.05 for the 
electric field and f{wt) = 2cos(2 wt) for the crossflow. The concentration throughout the system 
was initially uniform using a 0.021 wt% solution of BHb. The syringe, which served as the inlet 
reservoir, was filled with a volume of 8 mL of protein solution, while the outlet reservoir contained 
approximately 40mL of the same solution. 

For each of the experiments, the frequency of oscillation was chosen to be 0.068 ^ and the 
stroke volume to be 1.5 mL. For a = 0.25, an electric field amplitude of 9.795 ^ was required, 
and for a = 0.5, the amplitude was doubled. In both experiments, the concentration of the outlet 
reservoir was monitored with time. The results of the experiment using o = 0.5 is shown in Figure 
4. Note the concentration of the outlet reservoir increased by over 10% for 20 minutes of run time 
while the syringe concentration (not shown) was almost entirely depleted. The results of these 
experiments indicate that even with initially uniform concentrations, a net flux could be obtained 
for both choices of a, and furthermore, the magnitude of the velocity qualitatively agreed with 
theory. 

Experiments will continue with single protein species using various oscillatory shear flows and 
dimensionless mobilities under both transient and steady state conditions. From these experiments 
we hope to achieve quantitative comparisons of the net velocity and effective diffusivity of various 
protein species. Additionally, we will conduct multi-protein experiments in an attempt to verify 
that binary oscillatory crossflow electrophoresis as an effective binary separation technique. 
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Figure 4: Plot of dimensionless weight percent of Bovine Hemoglobin in the outlet reservoir versus 
time for a dimensionless electrophoretic mobility of on — 0.5. 
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